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Abstract 

We consider the quantum state describing the Disoriented Chiral Condensate (DCC), 
which may be produced in high energy collisions. We show how a mean field treatment 
of the quantum equations corresponding to the classical linear sigma model leads to 
a squeezed state description of the pions emerging from the DCC. We examine various 
squeezed and coherent state descriptions of those pions with particular attention to charge 
and number fluctuations. We also study the phenomenology of multiple DCC domains. 
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1 Introduction. 

The possibility of Disoriented Chiral Condensate (DCC) production in high energy 
hadronic or heavy ion collisions is a subject of much current activity. After the first 
papers |IJ-[|J this area has attracted much interest in the last 2 years and many different 



aspect have been discussed in numerous papers [p|-[|19fl- The aim of this paper is to discuss 
the possible quantum descriptions of the disoriented chiral condensate, paying particular 
attention to its decay into pions. 

It is well known that the QCD Lagrangian is invariant (only approximately if nonzero 
masses for the light Nj quarks are taken into account) under global chiral SU(Nf)£ x 
SU(Nf)n, where Nf is the number of the light flavors. This symmetry is spontaneously 
broken to vector SU(Nf) v which leads to Nj — 1 (quasi)goldstone bosons - pions (if 
Nf = 2) or pions, kaons and rj mesons (if Nf = 3). The order parameter for this 
breaking is the vacuum expectation value of the quark condensate < ipip >. In our 
normal vacuum, that condensate is an isotopic spin zero condensate (the a). However one 
can imagine that under some special conditions, for example after high-energy collision, 
there is a region in which the condensate points in some other, arbitrary direction in 
isospin space. This "cool" region of disoriented chiral condensate, (DCC), would be 
surrounded by a "hot" relatively thin expanding shell, which separates the internal region 
from ordinary space. This picture was suggested by Bjorken, Kowalski and Taylor |J 
and is called the "Baked Alaska" scenario. After hadronization the interior disoriented 
vacuum will collapse, decaying into pions. The interesting signature of this process could 
be the coherent production of either nearly all charged or nearly all neutral pions. Such 



remarkable charge fluctuations have been reported in the Centauro cosmic ray events |20 
These events may be a manifestation of the DCC. Much of the theoretical impetus for 
discussions of the DCC have been based on classical models of QCD and of the classical 
pion field. But the pions that emerge from the decay of the DCC are quanta of that pion 
field, and hence we need a quantum discussion of the field emerging from the DCC. The 
principle results of this paper concern the possible forms of that quantum description. 

There are two important quantum aspects of the DCC. The first involves the formation 
and evolution of the DCC and the second its decay into pions. The evolution of the DCC 
has been discussed in terms of a simple, classical, toy model, the linear sigma model. 
We will show that a quantum treatment of that model, in the mean field or Hartree 
approximation, suggest that the pions emerge from the decay of the DCC in a squeezed 



state. Other quantum states for the pions coming from DCC decay have been proposed. 
The coherent state is the first that comes to mind for quantizing a classical field. It has 
been objected that a coherent state leads to large charge fluctuations. We will discuss 
various coherent and squeezed states and show that a natural outgrowth of a real, classical 
field theory o the pions is a real cartesian (in isospin) coherent or squeezed state and that 
these states do not have large charge fluctuations. We also discuss briefly how a few 
independent domains of DCC would affect the phenomenology. 

In Section 2 we briefly review the classical linear sigma model treatment of the DCC. 
This helps to establish our notation and connect our work to the DCC literature. In 
Section 3 we show how mean field quantization of the sigma model leads to squeezed 
states for the pions. Section 4 is a discussion of the various coherent state and squeezed 
state pictures that can be used to describe DCC decay. We discuss treatments that 
can be found in the literature and some new forms, with particular emphasis on solving 
the problem of charge fluctuations. Section 5 briefly discusses the phenomenology of 
independent domains of DCC and Section 6 presents a summary and conclusions. A 
discussion of the neutral fraction probability for classical but not spherically symmetric 
isospin distributions is given in Appendix A. Some technical material on squeezed states 
and their relationship to coherent states is presented in Appendix B. 

2 Simple Classical Model of the DCC 

Instead of considering the real hadron world with all it complexities, most discussions of 
the DCC begin with a very simple, toy, classical model describing the chiral dynamics, the 
linear sigma-model. This expresses the dynamics in terms of a four component classical 
field, 4> a = (a, 7?), where a and n are isoscalar and isovector fields, respectively. The action 
is given by 



S = d x 



l -d,rd^ a -\{rr-v 2 ) 2 + Ha 



(i) 



where H ~ m q describes the small explicit chiral symmetry breaking due to the quark 
masses. The pion mass is m^ = H/ f n = X(f 2 — v 2 ), where /^ =< o >. The sigma 
meson mass is m\ = SXf 2 — Xv 2 « 2A/^. With m n = 135MeV,m cr = GOOMeV and 
U = 92.5MeV one has § A = 20 and v = 87 AMeV . 



The fields <p a = (er, 7?) parameterize the quark condensate 

7,- , TT 4 .^^ 1 cr + in m + \ , nS 

< V>Vi >~ U] = cr + in? =[ . _ . 2 

J ■> \ —m a — m J 

Let us note that if one fixes the "absolute value" of the quark condensate, i.e. the 
determinant of the matrix W, det U = a 2 + it 2 = f 2 one can get another representation, 
which is used usually for non-linear sigma-models: 

(17TT \ TT TTT TT 

— = f n cos — + iU — sin — (3) 

Jn ) Jit TT J n 

where n = \tt\. In case of small n « f n =< o > one can easily see the equivalence of 
these two parameterizations. 

In the usual vacuum one has < o >= f n , <7r>=0 and since a is an isoscalar there is 
an isoscalar condensate < ip l ipj >~ <$'■ only. However one can consider another configura- 
tion - < a >= f n cos 9 and < 7? >= f n n sin 9, here n is some unit vector in isospace. This 
condensate corresponds to the DCC, i.e. some classical pion field configuration that points 
in an arbitrary direction in isospin space, and which is metastable. Rajagopal and Wilczek 
have pointed out that the long wave length modes of the condensate grow exponentially 
fast after a quench ||. After some time the metastable disoriented condensate decays into 
the normal vacuum by emitting pions. All directions in isospin space are equally probable, 
thus there is reasonable probability that the decaying DCC will lead to an event with an 
unusualy large number of either neutral (71°) or charged (w^) pions. That would be a dra- 
matic experimental signal for the DCC, a signal that may already have been seen in the 
Centauro events |[20|| . Using the classical picture of DCC, which has equal probability for 
all isotopic orientations, one finds [l] - || the probability P(f) = l/(2y/f), Jq P(f)df = 1 
for the neutral pion fraction, / = N^o /{N pi + +N n - +N n o). Thus this picture leads to a far 
higher probability for finding nearly all charged pions, / ~ or nearly all neutral pions 
/ ~ 1 in the decay of the DCC, than the very small probability for such distributions 
that would be found in a purely statistical decay of an excited pionic region. 

The simple, classical, linear sigma model helps to motivate the DCC and its remarkable 
prediction for the distribution of pion charges in its decay. However the model is purely 
classical while the decay products of the DCC are pions, the quanta of the field. It is 
therefore important to study the quantum evolution of the DCC and its wavefunction 
for decay into pions. The usual approach to quantizing classical fields is to use coherent 
states pl| . We will discuss the various forms of coherent states that might be appropriate 



to the decay of the DCC and some of their problems, advantages and phenomenological 
consequences in Section 4, but before turning to a description we examine models of the 
quantum evolution of the DCC. We will see in the next Section that a quantum theory of 
the evolution of the DCC treated in the Hartree approximation leads to a picture of the 
decaying DCC in terms of two-pion coherent states, the so called "squeezed states." 

3 Quantum Evolution of the DCC and Squeezed States 

As has been shown in ]TJJ the most natural quantum description of the DCC is given in 
terms of squeezed states and the pion production in this state is a nice analog of an ampli- 
fication of oscillations in a quantum parametric oscillator. To get the squeezed quantum 
state for the DCC let us consider the mechanism for amplification of the long wavelength 
pion modes, suggested by Rajagopal and Wilczek in paper ||, where the the dynamics 
of the 0(4) linear sigma model after quenching was considered. The amplification of the 
long wavelength pion modes was found in the period immediately after quenching. This 
amplification leads to the coherent pion oscillations, i.e. to the creation of the DCC. Such 
a behavior can be understood if one considers the equation of motion for the pion field H 



|^(/c, t) + \P + A(< 2 > it) - v 2 )}7r(k, t) = (4) 

where we substituted (as in ||) the </> a a in the nonlinear term in (Q) by its spatial average 
< (f) 2 > (£). This is nothing but the Hartree or mean field approximation. If in the initial 
conditions one has < 2 > < v 2 then the long wavelength modes of the pion field with 
k 2 < X(v 2 — < (p 2 >) start growing exponentially. The < 2 > (t) starts to oscillate 
near the vacuum expectation value < a > and after some time the oscillations will be 
damped enough so that all the modes will be stable. Thus we see that at the classical 
level, each long wavelength mode is described by the equation for a parametrically excited 
oscillator and one gets the DCC as a result of the amplification of the zero-point quantum 
fluctuations of the pion field. 

This picture is similar to one which was discussed in [p2 |, where relic graviton produc- 
tion from zero-point quantum fluctuations during cosmological expansion was considered. 
For graviton mode with momentum n the equation y" + [n 2 — (R" /R)]y = was obtained, 
wherei?(^) is the scale factor of the metric ds 2 = R 2 (r])(dr] 2 — dx 2 ) and a prime represents 
d/drj. One can see that this equation is equivalent to the pion equation (f|) if the scale 
factor R is connected with < 2 > (t) as X(v 2 - < <fi 2 > (£)) = R" / R. 



Our problem now is to present the quantum mechanical formulation in terms of pion 
creation and annihilation operators and to get the wave function of the DCC. In the 
mean field approximation the wave function |\l/ >= T7j kl^ > ik * s ^ ne P r °duct of the wave 
functions for each mode with momentum k and isotopic spin index i. Later we shall omit 
%. The equation of motion (f|) for each mode 7v(k,t) can be derived from the Lagrangian 

L k = h 2 (k,t)-±n 2 (k,t)n 2 (k,t) 

n 2 (k,t) = P + \(< ( f ) 2 >(t)-v 2 ) (5) 

The wave function \ip >£ obeys the Schrodinger equation 

i^>% = H k (t)\1» % = [\vl + \n 2 {k^ 2 {k)]\i'>k (6) 

where 7r(k) and V^ = —id/dir(k) are the quantum- mechanical coordinate and momentum 
for the mode with the spatial momentum k. 

After finding the wave function \ip >£ and constructing the full Hartree function |\1/ > 
one can get the value of the mean field < <fi 2 > (t) =< \I/|o" 2 + 7r 2 |\l/ > and then obtain an 
equation for the mean field evolution. The last part of this selfconsistent approximation 
will not be discussed in this paper, our aim is to analyze the general features of a Hartree 
wave function, taking the mean field < <ft 2 > (t) as granted. 

One can rewrite the Hamiltonian in (|6|) in terms of creation and annihilation operators 
which make it diagonal at any given moment. It is evident that we are interested in getting 
the wave function in terms of creation and annihilation operators of ordinary pions, so we 
must diagonalize the Hamiltonian at t — > oo, when the oscillation of the < <fi 2 > (t) will 
be damped. Thus we define 

a (k) = ^Mk) + iv^ ot ^ = ujicMk)-^ (7) 



2u(k) V 2w ( fc ) 



where u(k) = Q(k, oo) = \Jk 2 + m 2 . It is easy to see that the Hamiltonian (up to a 
constant) is 

H k = Uk)[l + ^£-W(k)a(k) - " 2 W-^ M V (fc) + * j2 m (8) 

It is the a 2 (k) and a* 2 (k) terms in the Hamiltonian that transform the initial vacuum |0 > 
into a squeezed state S(r, 0)|O >. The squeezed state is given by (pi) 



K>=5(O|0>=exp 



\ (tf - e(^) 2 



|0 > (9) 



and we are writing £ = rexp(i0). More generally, one can consider the squeezed coherent 
state S(£) exp(aa^)|0 >. Some technical details about squeezed states are presented in 
Appendix B. 

To calculate the squeezing and phase parameters r and let us consider the solution 
of the Schrodinger equation (^j) in the coordinate (here it is 7r) representation ( we omit 
the mode label k for a moment) 

< vr|^ > (£) = C(t) exp(-B(t)n 2 + D(t)w) (10) 

where we also take into account the possiblity of the time-dependent expectation value 
of the 7r field < n(t) >=< ip(t)\Tc\ijj(t) >~ D(t)/2B(t). This state corresponds to the 
squeezed coherent state. For B = uj/2 this wave function describes the usual coherent 
state and, if D — 0, it is the vacuum state. For all other values this wave function 
describes the squeezed state ([J) where the parameters r and are connected to B by the 
relation |4f (see also || - ||) 

u coshr + exp(2z0) sinhr 



B 



2 coshr — exp(2i0) sinhr 



uj 2 + A\B\ 2 . nl 1 ImB 

cosh2r = ; sin 20= — (11) 

AuoReB ' Y sinh2r ReB v ; 

Substituting (|10| ) into @ one gets equations for B(t), C(t) and -D(t): 

•J3(t) = 2B\t) - lfi 2 (t), i£>(«) = 2B(t)D(t), zgj| = fl(t) - ^ 2 (i) (12) 
which means that B(t) is related to the solution of the classical equation (|J) 

m = ^m> m+tf(t)m = o as) 

in terms of which D(t) can also be obtained 



D(t) = £>(0)exp(-2z^ B(r)dr) = D{0) exp f - jf ^rj = 0(0)^ (14) 

and C(£) = C(0)exp (if (B(r) — D 2 (r)/2)dTj is an evident phase factor. The equation 
fljjD can be viewed as a Schrodinger equation describing the wave function ip(t) of a 
"particle" with mass m = 1/2 on a line t having energy u 2 (k) = k 2 + m 2 and moving 
through the potential barrier V(t) = — A(< 2 > (£) — / 2 ) 

- ^ + A(/ 2 - < 2 > (t)ty(t) = u 2 {k)i,{t) (15) 



*- 


-»■ +00 








t- 


-> —00 






(16) 


al 


1 c+12 


1 c+12 

\ D R\ — 


1 O — 12 


1 C-|2 



Far from the barrier, i.e. at t — > ±00 one has V(±oo) = and the general solution of the 
Schrodinger equation at t — > ±00 is a superposition of left and right moving waves 

^+(t) = S+e~ iuj{k)t + S+e +luj(k)t , 
</>(£) = S R e- w(k)t + SZe +iuj{k)t , 

Due to unitarity the total fluxes at t — ■> ±00 must be equal (S^ 
and one can find 

Sjj = cosh rSft — e 2ie sinh rSj^ 

S^ = -e- 2ie sinhrSx + coshrSl (17) 

where 9 is the scattering phase and the factor r is defined by the probability of transition 
through the barrier. 

Let us remember that if we are starting from the vacuum at t — > —00, i.e. from 
B = u/2, one must have Sr = 0. Even if we assume that the initial state is not the 
vacuum, but some coherent state with nonzero -D(O) we still have to put B = u/2 initially. 
This means that at the left (large negative t) we have only a left moving outgoing wave 
S'^e +4w ^- )t . At the right (large positive t) one has both left and right moving waves, i.e. 
the incoming g+g+^w* anc } reflected S^e - " 1 ^* waves. The transition coefficient can be 
obtained from (|X7|) by puting S^ = 

! -' - (18) 



St I 2 cosh r 



'L 



Now let us calculate B(t) = —(i/2)(ip/if)) at large positive t. Using QT5| ) one can find 
after simple calculation (restoring the k dependence) : 

u(k) cosh r(k) + exp[2i(9 — u(k)t)] sinh r(k) 

^ ' ~ 2 coshr(k) - exp[2z(^ - u(k)t)] sinhr(fc) ^ ' 



in complete agreement with (pT), with the phase factor <j>{f) depending on time as = 

e-u(k)t. 

The squeezing parameter r(k) depends on the absolute value, k, of the mode spatial 
momentum k and is determined by the probability of tunneling through the potential 
barrier V(t) = A[/ 2 — < 2 > (£)]. Let us note that tunneling takes place precisely when 
uj 2 (k)— V(t) < 0, i.e. when the classical long wavelength modes are exponentially amplified 
and we see once more that squeezing is ultimately connected with the exponential growth 

8 



of the classical long wavelength modes and the squeezing for each mode k is determined 
by the function < 2 > (t) . This is the only input information we must know to calculate 
the DCC wave function. 

After calculating B(t) one can study the effect of a displacement D(t). All we need 
to know is the ratio ip(t)/ip(0), where "0" is an arbitrary initial time and we shall put it 
to t- — ► — oo. It is easy to see that (we put the scattering phase 9 = 0) 

S£ exp(+iu(k)t) 



D(t+ -»• oo) = D(t- -»• -oo) 



Sr exp(—iu(k)t+) + S£ exp(+iu(k)t^ 



e -iu>(k)(t+-t-) t (20) 

coshr (1 — exp(— 2iu(k)t+) tanhr) 

and one can see that at large t+, on average, D(t+) is small due to the suppression factor 
1/ coshr. Only at moments t = im/uo during the short intervals St ~ e~ 2r /u> does one 
get large D(t). However these pulses may be an artefact of the approximation, with the 
average D(t) being of order e~ r . Thus we see that for those modes for which the squeezing 
is large, the initial displacement, i.e. the coherent part of the wave function, disappears 
at late times independently of how large it was in the beginning. The mean value of the 
pion field < ir(t) >~ D(t)/B(t) and one gets 

< ir(t) >~ e -M fe X*+-*-) — J: ^_ (21) 

coshr (1 + exp(— 2iu{k)t+) tanhr) 

and on average, < n(t) > is also suppressed by a factor 1/ coshr. Only at moments 
t = (nn + 7i/2)/uj during the short intervals St ~ e~ 2r /u does one get large < 7i(t) >. 
However these pulses may be also an artefact of the approximation we use as the pulses 
for D(t + ) and on average the classical pion field is small. This is a very important point 
demonstrating that for low momentum pion modes, the classical pion field, even if it 
were produced at some intermediate step, will ultimately be negligible. This is for the 
displacement field. The average of the square of the pion field, < (7r(t)) 2 > is, of course, 
not small but is the exponentially growing condensate of the DCC. This is the part that 
leads to the squeezed state. 

In the quasiclassical approximation it is easy to calculate the squeezing parameter 
r(k): 



r{k) = 2Re / dt^\[v 2 - < 2 > (£)] - k 2 (22) 



which is valid for small k when r(k) » 1 The average number of particles in each mode 
is < Nk >= sinh r{k) and we see that < Nk > sharply decreases with the increase of k. 

9 



To make a rough estimate let us consider a simple model for < <ft 2 > (£). We assume that 
< <p 2 > (t) = in the interval t e (0, r) and is equal to its usual v.e.v. < 2 > (£) = f 2 
outside this interval. Such behavior for < 2 > (t) is a very rough approximation to a 
more realistic behavior obtained in |5| in numerical experiments, nevertheless one can use 
it for some preliminary estimates. The transmission coefficient in this case is well known 
(see for example |28[]) and in our notation takes the form 

1 _ 4u 2 (k)(ml/2-u\k)) 



cosh r(k) Auj 2 (k)(ml/2 - uj 2 (k)) + (ml/2) 2 sinh 2 rJml/2 - u 2 (k) 



where we used the relation m 2 = 2Xf 2 . Then the average number of particles with 
momentum k is 

/ m 2 \ 2 sinh 2 t J ml/2 - u 2 (k) 
<N(k)>= sinh* r(k)=[-f) ^ (k) J m { /2 _^ {k)) (24) 

One can estimate the average number of particles < N(k) > for different k. For 
r we shall use an estimate r = 3 — 6 m^ 1 which comes from the results of || r = 
1 - 2 (200MeV)-\ Then for k = one gets < JV(0) >^ 10 2 - 10 3 , at k = m^ the 
< N(m n ) > is approximately twice as small but at k — 3m n ~ 400MeV it is two 
orders of magnitude smaller < N(3m n ) >m 1 — 10. Of course this is the very crude 
estimate, however it demonstrates the qualitative features of the phenomenon - the sharp 
exponential dependence on k and large amplification factor in r which is of order of mpT, 
where r is the characteristic time of damping of < 2 > (t) oscillations. 

In modeling a realistic "Baked Alaska" scenario one should include the effect of the 
expansion. In ||] it was suggested that the expansion be described by including the 
term air/a in the equations of motion, where a(t) is a scale factor for the expanding 
plasma. This is the same as considering the problem in a space-time with metric ds 2 = 
dt 2 — a 2 (t)dx 2 . It is easy to show that choosing the new time t(t) = J dt/a 3 (t) one gets the 
same Schrodinger equation as in (^) but with new Q 2 (k,t) = a 6 (t)Q 2 (k,t). Thus taking 
into account the expansion one gets squeezed state again but with parameters which 
depend on the details of an expansion - the scale factor a(t). In a case of anisotropic 1 + 1 
expansion, which seems more realistic at early stages of the DCC formation (if any) and 
which was considered in [0] , || and later in |1| the convenient coordinates parameterizing 



the interval are proper time r and rapidity rj 



1 , ft — x 



V¥^, V=-^{—) (25) 
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and ds 2 = dr 2 — T 2 drf — dx\ In this case one gets the Schrodinger equation for evolution 
in r with the fi 2 (k±, k v , r) where momentum k v appears in the combination k 2 /r 2 , which 
means that at small r one has amplification only at very low k v . One can get analogs of 
(|2~3|). In this case it will be N(k v , k±) and one can get the distribution in rapidities r\ after 
taking the Fourier transform in k n . 

Thus we see that quantizing the linear sigma model for the DCC and treating the 
non-linear equation in the Hartree or mean field approximation leads to a simple picture 
of the quantum DCC that suggest a squeezed state treatment of the emerging pion waves. 
We also find the exponential amplification of the long wavelength modes that is seen 
classically leading, in a simple model, to a very large number of low energy coherent 
pions. We now turn to various treatments for the wave function of those emerging pions. 

4 Quantum Wavefunctions for Decay of the DCC 

There are many ways to go from a classical pion field theory to a quantum state, and 
many of these have appeared in the literature of DCC and other places. In this Section, 
rather that only give the squeezed state description we developed in Section 3, we present 
and contrast these various treatments. Thus this section does not depend on the squeezes 
state scenario, but rather discusses a range of possible quantum treatments of the classical 
decaying DCC. 

The simplest quantum description of a classical field is given in terms of a coherent 
state |H[ 



|tt >= exp [ ) exp ( f d 3 kn(k) ■ sJ(k)\ |0 > (26) 

where n is the isotriplet pion field and and at are corresponding creation operators. The 
mean number of pions in the state is given by 

N = f d 3 kir ■ it* (27) 

which number also comes into the normalization of the coherent state. The root mean 
square number fluctuation in this state is V N, and hence the fractional number fluctuation 
decreases for large N. If one defines position and momentum like variables for each mode, 
k, by 

f= 7! (at + a) (28) 

11 



and 

P=-^(a f -a) (29) 

one also finds that the coherent state is a minimal uncertainty state. It is these two 
properties, minimal uncertainty and small number fluctuation in the large N limit, that 
make the coherent state attractive as the quantum representative of the classical field. 

However this coherent state wave function has been criticized in Q because, for arbi- 
trary ir(k) such a description leads to large charge fluctuations in the pion system. In the 
remainder of this section we will discuss ways of addressing this criticism while maintain- 
ing as much as possible of the advantages of the coherent state, minimum uncertainty and 
small number fluctuations, and also maintaining the relatively large probability of finding 
a neutral pion fraction near zero or one. We will divide our discussion in two parts, states 
with identically zero charge and charge fluctuations, and states with only average charge 
zero. 

Let us begin with states with average charge zero but still having charge fluctuations. 
Consider the coherent state defined in fl2"E|). Let us use a cartesian rather than a spherical 
tensor representation for the isospin of the pions, both for the creation and annihilation 
operators and for the classical field, ir(k). Since the classical field equations are real, we 
might expect that classical field configuration to be real as well. Hence a real cartesian 
quantum description is a natural outgrowth of the classical beginnings of the DCC. We 
will now show that this leads to zero expected charge in the coherent state and small 
fractional charge fluctuations. In terms of cartesian pions, the number operator for tt + is 

given by 

1 , i 

n + = ~(a\a x + a^cLy) + -{a y a x - a x a y ) (30) 

and for 7r~ by 

1 i 

so that the charge, Q = n + — n_ is 

Q = i(a y a x — a x a y ) (32) 

This is for each mode, k, the total charge is the integral over k. It is easy to see that, for 
real cartesian pion field, the expected value of the charge in the coherent state is zero. 
The expected value of the charge squared is 

< Q 2 >= J d^l + „J) (33) 
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which is a number of order N. Hence the fractional root mean square charge fluctuations 
are of order 1/viV and therefore small in the large N limit. It is clear that the real 
cartesian coherent state also continues to be a minimal uncertainty state. 

We would expect that in the DCC the classical pion field will point in some fixed 
direction in isospin space for all k. Suppose we label that direction with the usual polar 
angles, 9, <p, and continue to use the real cartesian coherent state. Then it is easy to see 
that the average number of 7r°'s is given by A cos 2 6*. The neutral pion fraction is thus 
given by / = cos 2 #, and corresponds to a probability of P(f) = 1/(2a/J), as before. The 
root mean square pion charge fluctuation is given by viVsin^ and the fractional charge 
fluctuation by sin9/vN, which is small in the large N limit. Hence the real cartesian 
coherent state has many desirable features for describing the radiation of the DCC. Its 
major faults are that it does not have identically zero charge and that it is not well linked 
to a simple quantum theory of the DCC dynamics itself, which suggests squeezed states. 
|19| We now turn to these issues. 

Let us begin with the application of squeezed states to the problem of decay of the 
DCC. We first review some features of squeezed states. For a single mode created by the 
operator b\ the normalized squeezed state may be written 

\a >= (1 - |a| 2 ) 1/4 exp (f (& f ) 2 ) |0 > (34) 

where a is a complex number with modulus less than one and is related to the parameters 
of (Bp. (See Appendix B for more details.) The squeezed state is a minimum uncertainty 



state, for position and momenta defined in terms of the operators b as in ( p8|) and (29). 
In the case of the coherent state, that minimum uncertainty is achieved by having each 
of Ax and Ap take on their minimum values. For the squeezed state, only their product 
need be minimum, and one can be large while the other is "squeezed" to keep the product 
fixed. An important difference between the squeezed state and the coherent state comes 
in the number fluctuations. For the squeezed state we find 

<a|(6 t 6) 2 -(A) 2 |a>=2A(A + l) (35) 

which means that the fractional root mean square fluctuations are of order 1, in the large 
N limit rather than of order 1/viV as they are in the coherent state case. 

In spite of the large number fluctuations in the squeezed state, let us look at its 
application to the DCC. Suppose we construct the operator tf from a real isospin rotation 
of the cartesian pion operators 



b* = cos 9a\ + sin 9 cos <\>a\ + sin 9 sin (f>a) y (36) 

13 



where 9 and are the usual polar angles. We may now label our squeezed state \a, 9, (f> >. 

_ I |2 

It is clear that the mean number of pions in the state is still iV = 1 _jL 2 - The expectation 
of the total charge in this state is easily seen to be zero. The expected value of 7r°'s is 
given by 

< a,9,(p\ala z \a,9,4>>= N (cos 9) 2 (37) 

The neutral pion fraction is / = cos 2 9 and its probability is again P(f) = — j=- Hence 
the squeezed state with arbitrary real cartesian isospin direction has zero average charge, 
minimum uncertainty, and a large probability for neutral fraction near zero or near one. 
However it has large (of order N) charge fluctuations. 

To deal with the charge fluctuations, it was suggested in M that the quantum state 
that describes the DCC be an isosinglet. A particular isospin zero wave function was 
considered a long time ago by Horn and Silver [23[. Consider the operators A and A\ 
with 

^ = 2aW- ~ (4) 2 = "((4) 2 + (al) 2 + (at) 2 ) (38) 

which is clearly an isoscalar operator. Any function of A* operating on the vacuum will 
create a state with isospin zero and thus an eigenstate of charge with zero eigenvalue. 
Such a state will have not just zero average charge, but zero charge fluctuations as well. 
We note, in passing, that A, A^ and N = J2i a \ai are the generators of an sl(2) algebra 

[A, A^} =AN + 6, [N, A] = -2A, [N, A ] ] = 2A f (39) 

A particularly simple state of isospin zero and exactly 2N pions has been suggested 
by a number of authors, cf. j||. It is of the form 

|tf >= ; l (A) N \Q > (40) 

y/(2N+l)\ 



Using Stirling's formula for the asymptotics of n\ ~ V2raexp(nln(n/e)) it is easy to see 
that the probability of having 2n neutral pions in the state |\l/ > is |p3|1 , Q 

f/Vn 2 2 27V (2nV i 

P(n, N) = -Kr T7 , , „ ~ \/N/n, n,N»l (41) 

v ' ; (2N + 1)\ (n\2 n ) 2 V / ' ' v ) 

and corresponds to the 1/a/J distribution of the classical picture. 



The same distribution would be obtained if an arbitrary relative phase factor were 
inserted between the charged and neutral creation operators ( 2a' + a_ — exp(i9)(a' 



2> N 
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This might seem to imply that the zero isospin condition is not so important, and that 
one can introduce any arbitrary factor. This is not the case. With any arbitrary factor in 
front of (aj) 2 , the states created will still have the third component of isospin equal to zero 
and hence be zero eigenstates of charge, but the operator A* will now be a combination of 
isoscalar and second rank iso-tensor All directions in iso-space will not be equivalent and 
we will no longer have a 1/y/J neutral fraction distribution. For example if the arbitrary 
factor is zero, there will be no neutral fraction. 

The state |\P > considered above is quite restrictive. It has exactly 2N pions. A 
quantum state arising from a classical field is not expected to have a sharp number of 
quanta. We now consider two isoscalar states that are closer in spirit to the coherent 
state, and have an indefinite number of pions. First consider the squeezed state made by 
exponentiating A'. We write 

\q,a>=(l-\a\ 2 )iexp(^A^\0> (42) 

To understand the normalization factor note that this state can be written in terms of 
cartesian pions as a product of three simple squeezed states, 

|*, a >= (1 - H 2 )! exp (^((4) 2 + K) 2 + (4) 2 )) |0 > (43) 

The probability of finding 2n neutral pions in this state is 

OO 1 1 

P(n)=^ JW P(n,N)^ 7s (44) 

- again the same distribution. Although this state is a minimum uncertainty state, and 
there are no charge fluctuations in it, there are large particle number fluctuations. The 
number operator can be written in terms of the cartesian number operators as n x +n y +n z . 
Its expectation value in the state \^f,a > we call N. Using the product nature of the 
state we can then calculate 

a 2 =< ^,a\(n x + n y + n z ) 2 \^,a > -N 2 (45) 

using the results for a simple squeezed state. We find 

N 2 
a 2 = — + 2N (46) 

This shows that the number fluctuations grow like N, or the fractional fluctuations like 1 
rather than like 1/viV. 
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As an alternative for the iso-scalar state one can consider the iso-spin projected co- 
herent state first introduced by Botke, Scalapino and Sugar, |30| and recently used to 
describe the coherent pions emerging from nucleon-antinucleon annihilation. |5l| , f32|j 
The 1 = projected coherent state is 

| A, >= AT J dfe md3ka l- f \0 > Y* (f) (47) 

where the normalization is given by 

M={A-n ]G {-iN)y l l 2 (48) 

and the mean number of particles, N is given by 

N = [ d 3 k\(k)X*(k) (49) 



This is an iso-scalar state and hence has zero average charge and zero charge fluctuation. 

It is also easy to show that for large N, the variance, a grows like vN so that the 

fractional particle number fluctuations decrease for large N. However the state |A,0 >, 

is not a minimum uncertainty state. To study the uncertainty in the state let us look at 

the uncertainty in the third component of isospin (all components are equivalent since we 

are in an / = state). Define 

i 

£o = -y=(a + al) (50) 

and 

Po = -i-T=(a Q -al) (51) 

where the subscript zero refers to the third component of isospin. It is easy to see that 

< A,0|x |A,0>=< A,0|j9 |A,0 >= (52) 

This is also a feature of the 1 = squeezed state. We need to look at x 2 and p 2 . They 
are given by 

4 = -{a 2 + {al) 2 ) + ala + l/2 (53) 

and 

pl = --(a 2 + (al) 2 )+ala + l/2 (54) 

One finds (taking A(fc) as real) 



< A, 0|oo|A, >=< A, 0|(aJ) 2 |A, >= — (55) 
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and 

<A,0|ata |A,0>=pll^ (56 ) 

3 jo{-iN) 

Finally in the large N limit one obtains 

TV 
AxAp = — (57) 

Hence the isospin projected coherent state is not a minimal uncertainty state. 

Thus we see that there are many classes of trial wave functions that can reproduce 
some of the phenomenologically desirable features of DCC decay. It is natural to ask the 
question what is the most natural class of these functions and what are the dynamical 
mechanisms leading to the generation of these functions. In the previous section we ar- 
gued that simple quantum dynamics for the evolution of the DCC suggests the squeezed 
state as being most natural. In this section we have presented a number of quantum state 
candidates for discussing the pion state from the decay of the DCC. We have examined 
them from the point of view of charge fluctuations, number fluctuations, minimal uncer- 
tainty and neutral fraction probability. None of the candidates pass all criteria, but the 
cartesian pion squeezed state and the cartesian pion coherent state seem the most attrac- 
tive. The first has the better dynamic credentials, but it also has much larger number 
fluctuations than the second. 



5 Domains of DCC 

We have seen that most reasonable quantum treatments of the pions from decay of the 
DCC lead to a probability P(f) of neutral pion fraction / of P(f) = 1/(2-^//). This 
means that the probability of finding nearly all charged pions or nearly all neutral pions 
is quite large compared with what it would be for a purely random collection of N (JV 
large) pions. The probability of finding / within e of zero (corresponding to nearly all 
charged pions) is y/e, while the probability of finding / within e of one (nearly all neutral 
pions) is \e for small e. It is these surprisingly large probabilities that has fueled much of 
the interest in DCC, since they imply a strong and unambiguous experimental signal for 
the DCC. But in reactions in which a DCC can be created it is possible that not a single 
domain of DCC will be created but that perhaps a number of disconnected domains will 
be created each with a separate random direction in isospin space for the condensate. 
What then will be the probability of an overall neutral pion fraction /? 

17 



Suppose there are a number of regions, domains, or bubbles of DCC and that the 
isospin condensate direction in each is arbitrary. If, of course, there are many many 
such domains the average of / over the domains will be 1/3 with very small variance. 
In particular the chance of finding / near zero would be very very small. But a more 
likely scenario is that there may be a few domains, more perhaps than one, but not many 
many. Then the probability of finding small / is not prohibitively small. For example, 
suppose there are two domains, and we average the pion observations over both. Then, 
the probability of finding neutral fraction / is it/2 for / < 1/2 and 7r/2 — 2 arccos( 



for 1/2 < / < 1. Similar forms can be easily obtained for n domains. In particular one 
sees that the probability of finding neutral fraction less that e is y/e for the one domain 
case, and (7r/2)e in the two domain case. For n domains the probability of finding neutral 
fraction less than e goes like e?. For reasonably small numbers of domains this is not 
ridiculously small. 

Hence one might build a phenomenological model of how the DCC manifests itself in 
terms of a small number of statistically independent domains of DCC. That might be 
closer to what happens in heavy ion collisions, say, and yet holds out hope of showing an 
interesting signal. 



6 Summary and Conclusions 

The idea that a Disordered Chiral Condensate (DCC) might occur in a rapid quench after 
a high energy collision is an attrative one. The signature of a DCC would be the emission 
of a large number of pions either nearly all neutral or nearly all charged. Such a remarkable 
process may have already been seen in the Centauro events. [ED| Most discussions of the 



DCC concentrate on how it arises in a simple classical field theory, the linear sigma model, 
and how the long wave length modes grow until they encounter the hot boundary of the 
quenched region ("Baked Alaska"). But the particles detected from the DCC are pions, 
the quanta of the field, and hence a quantum discussion both of the growth of the DCC 
and of its deacy is needed. In this paper we have attempted to provide the beginnings 
of such a discussion. By quantizing the linear sigma model in the mean field or Hartree 
approximation, we have shown how the quantum DCC grows and how it leads to a 
squeezed state of the pions. More generally we have examined quantum wave functions 
for the pions from DCC decay, both of the coherent state and squeezed type. We have 
shown that the classical starting point leads naturally to a quantum description of the 
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pions in terms of real cartesian (in isospin) pions, and that that description avoids many of 
the problem of charge fluctuations in the final pion state. We have also tried to emphasize 
the phenomenological differences among the various quantum descriptions. We have also 
discussed how multiple independent domains of DCC, might manifest themselves. These 
are likely to occur in heavy ion collisions. 

Much remains to be done in this growing subject. Better classical and corresponding 
quantum models are needed, with a better understanding of how the classical models 
connect to particular quantum descriptions. The notion of solving the difficult strong 
interaction dynamics classically and then tying the solution onto a quantum state is a 
very promising one and has already shown some success in the description of nucleon 
antinucleon annihilation. [J31J1 , |32| However the real test of the DCC concept must come 
from experiment, and therefore more phenomenological as well as purely theoretical work 
is need as we approach the time of data. 

The work of R. D. Amado is supported in part by the U. S. National Science Founda- 
tion. The work of I. Kogan was supported U.S. National Science Foundation grant NSF 
PHY90-21984. 

Appendix A: Ellipsoidal Isospin Distributions 

All discussions of the neutral fraction proceed from the assumption of a spherically sym- 
metric a priori distribution in isospin for the DCC condensate. For completeness we 
investigate in this appendix what happens if we replace this with an ellipsoidal distribu- 
tion. We continue to assume that the neutral pion fraction, /, is given by cos 2 9, but now 
ask for the probability of finding a given / over an ellipse rather than a sphere. To do this 
take ordinary cartesian coordinates in isospin space. Then cos 2 9 = x2 2 2 2 . We must 
take an ellipse that is symmetric in x and y to preserve zero charge, hence we take for 

our ellipse 

x 2 + y 2 z 2 



a 



2 



+ pr = l ( 58 ) 



The probability of neutral fraction /, P(f) is then 



Pif) = N j K1 -t+£- £ w _ _*_ )dxiyiz m 



x z + y z \ !fJ . z A 

x 2 + y 2 + z l 
where the normalization is given by 



x 2 + y 2 z 2 



l/N =5(1 ^- - —)dxdydz (60) 



a 



b 2 



19 



It is clear that the probability is normalized, J P(f)df = 1. Since / and P(f) are dimen 
sionless, the result for P(f) can only depend on p = a/b. We find 

1 



P(f) 



P 



(61) 



2v7(a + /(p»-l))§ 

This is seen to reduce to the usual spherical answer for p — 1. It is also easy to see that 
for p near 1, it is quite close to the usual answer. For p very large, only very small / 
dominates (nearly all charged) and for p very small, / near one dominates (nearly all 
neutral). 



Appendix B: Squeezed and Coherent States 

Squeezed states have been known for a some time in quantum optics and measurement 
theory ( for a review of squeezed states see, for example |23|]~ WJ\)- ^ ne srm plest one- 
mode squeezed state is parametrized by the two parameters r and and can be obtained 
by acting on the vacuum with the unitary squeezing operator S(£) 

K >= 5(0 10 >= exp [~ (^ 2 - £(a f ) 2 ))] (62) 

where £ = rexp(i0) is the squeezing parameter. The mean number of quanta in the 
squeezed state is N = sinh 2 r. 



To see this let us note that using the squeezing operator S one can make the Bogoluy- 
bov transformation b = SaS^ , tf = Sa^S^ and after some algebra one gets: 



b = cosh ra + exp(i(p) sinh ra\ tf = exp(— i<p) sinh ra + cosh rev 

a = cosh rb — exp(z0) sinh rb\ a) = — exp (—i(j>) sinh ra + cosh raJ 

The new operator b is the annihilation operator for the squeezed state 

6|f >= bS\0 >= SaS^S\0 >= Sa\0 >= 

Then it is easy to see that 

N =< £|a f a|£ >= sinh 2 r < £|66 + |£ >= sinh 2 r 



Let us note that one can rewrite [|33 
exp 



5 (e« 2 - a^) 2 ) = exp («(0a f 2 ) exp (/5(Oa 2 ) exp ( 7 (£)( fl t 



[ a+^ 



(63) 



(64) 



(65) 



(66) 
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because K + = a' f a^/2, K_ = aa/2 and Kq = l/2(a + a + 1/2) are the generators of the 
SU(1, 1) algebra: 

[K_, K+] = 2K , [K , K ± ] = ±K± (67) 

As a result, the squeezed state can be rewritten as 



|f >=exp( 7 /2)exp [a{£)a! 2 |0 > 



(68) 



To find a and 7 one has to calculate the A dependence of the overlap between a squeezed 
and a coherent state < A|f >=< 0|e Aa |f >. Then 

/-* 
2 ' - v ■ - ■ / ' ■ 2 ' " " ' " rfA 2 



exp (7/2) =< 0|e >, a = \< 0\a 2 exp (a(0<J 2 ) |0 >= ]- exp (-7/2) - K ^ > \ X=0 (Q9) 



Using ( |63"] ) twice, and the fact that b\£ >= 0, < 0\a^ = and [exp(Aa),a^] = Aexp(Aa) 
one gets after straightforward calculations 



d< Ajg > 



-Ae^tanhr < A|f > 



(70) 



The solution of this differential equation is 



< A|£ >= exp f e i0 tanhr ) < 0|f > 



(71) 



To find < 0|f > one has to consider the derivative d < 0|f > jdr and again using (| 
twice and 6|f >= 0, < O^ = 0. We find 



d .„ d 

— <0f >= — <0 exp 
ar ar 



(exp(— i(j))a 2 — exp(i0)(cr ) 2 



\0>-- 



11 1 

-exp(— i(j)) < 0|a 2 |f >= -exp(+z^) sinh 2 r < 0|6 t 2 |f > — coshrsinhr < 0|f >= (72) 

sinh r— < 0|f > H — coshr sinhr(sinh r — 1) < 0|f > 



dr 



and finally 



— < 0|f >= — tanhr < 0|f > 
dr 2 



(73) 



The solution of this differential equation using the boundary condition, < 0|0 >= 1 is 

1 



<0\£>- 



(74) 



vcoshr 

The overlap between the squeezed and coherent state is useful for expressing the squeezed 
state as a superposition of coherent states, recalling that the coherent states are complete. 
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After a simple calculation we get 

a = — tanhr exp(i0), exp(7) = — - — (75) 

cosh r 

Introducing new variable a = — tanhr exp(i</>) one can rewrite the normalized squeezed 
state as 

\a >= (1 - |a| 2 ) 1/4 exp (-(a f ) 2 ) |0 > (76) 

where a is a complex number with modulus less than one. The mean number of quanta 
in this state, N is given by 

I |2 

< alalia >= ; — rr (77) 

1 — \a\ z 

The squeezed state is a minimum uncertainty state, for position and momenta defined 
in terms of the operators a: AxAp = 1/2, where a(aJ) = (x =F ip)/v2- In the case of the 
coherent state, that minimum uncertainty is achieved by having each of Ax and Ap take on 
their minimum values. Coherent states have minimal quantum noise Ax = Ap = 1/ 'v2. 
For the squeezed state, only their product need be minimum, and one can be large while 
the other is "squeezed" to keep the product fixed. For example for (f> = one has 
Ax = exp(— r), Ap = exp(r). 
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